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Quantized spin models 
- Impurity and defect levels 
■ Exchange interactions 



Abstract. - We present a large- A 7 variational approach to describe the magnetism of insulating doped semi- 
conductors based on a disorder-generalization of the resonating- valence-bond theory for quantum antiferro- 
magnets. This method captures all the qualitative and even quantitative predictions of the strong-disorder 
renormalization group approach over the entire experimentally relevant temperature range. Finally, by map- 
ping the problem on a hard-sphere fluid, we could provide an essentially exact analytic solution without any 
adjustable parameters. 
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The metal-insulator transition (MIT) in doped semiconduc- 
tors (DS) [1] is one of the most fundamental, yet theoretically 
less understood problems in condensed matter physics. Even 
aside from their pivotal technological role, the DS have long 
been recognized as a bellwether system for the study of quan- 
tum criticality at the MIT. Careful transport experiments have 
revealed sharply defined critical behavior, although with expo- 
nents inconsistent with early theoretical predictions [2]. 
■ What are the basic physical processes that drive this transi- 
tion and localize the electrons? Important clues have been pro- 
vided by the thermodynamic response on the insulating side. 
Here, no magnetic ordering has been experimentally observed 
down to the lowest temperatures, while both the spin suscepti- 
bility and the specific heat display signatures of randomly inter- 
acting localized magnetic moments [1,3]. This puzzling behav- 
ior was largely explained by the Bhatt-Lee (BL) theory [4] of 
random singlet (RS) formation, using a strong-disorder renor- 
malization group (SDRG) approach [5]. 

The remarkable success of the BL theory provides strong 
support to the early ideas of Mott [6], who first emphasized 
that strong Coulomb repulsion may localize the electrons by 
converting them into localized magnetic moments. According 
to this picture, the MIT in DS should be viewed as a disor- 
dered version of the Mott transition, a phenomenon dominated 
by strong correlation effects. An appropriate theory should 
then be able to describe both the local moment magnetism in 
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the insulator and the transmutation of these local moments into 
conduction electrons on the metallic side of the MIT. Unfortu- 
nately, the SDRG approach of BL, which was so successful in 
the insulator, is difficult to extend across the transition. 

The essential challenge, therefore, is to develop an alterna- 
tive approach to Mott localization in a strongly disordered situ- 
ation, one that at the very least can reproduce the RS physics of 
Bhatt and Lee. An attractive avenue to describe strong correla- 
tions has emerged in the last twenty years from studies of var- 
ious Mott systems, based on resonating-valence bond (RVB) 
ideas of Anderson [7] and others. At the mean field level, 
these theories provide variational wavefunctions for quasipar- 
ticle states, which become exact in appropriate large- N lim- 
its [8]. Very recent work has extended similar variational stud- 
ies to disordered systems, providing a description of phenom- 
ena such as disorder-induced non-Fermi liquid behavior [9], but 
did not address the physics of inter-site spin correlations central 
to the BL paradigm. 

In this Letter we examine an appropriate t-J model capable 
of describing the Mott transition in a disordered environment. 
While the large- N limit of this model generally reduces to an 
RVB-like variational problem, here we concentrate on the lo- 
calized (i — ► 0) limit in the presence of strong positional disor- 
der modeling the insulating DS. We show that: (i) the large- 
formulation quantitatively reproduces all the key features of the 
RS regime; (ii) an accurate analytic solution of the variational 
problem can be thus obtained, providing closed form expres- 
sions for various physical quantities; and (iii) the approach can 
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be directly extended to the metallic side, eliminating the main 
stumbling block in attacking the MIT in DS. 

Model and large-N formulation. We start with the large-A 
formulation of the two-orbital t-J model, 

H = ]>> k + e )cl a c ka + £ U 3 flf ja (1) 

+^ E ■ s i + -7^ E ; '' kr c L/- + H - c -), 

under the constraint of no double occupancy on the /-orbital. 
Here each lattice site corresponds to a donor or acceptor which 
is randomly distributed in a periodic -boundary 3D cube of vol- 
ume Vq = Nq/pq, where Ao is the number of dopant sites and 
po is the doping concentration. We stay at half-filling for the 
uncompensated DS, ^ k CT c L c k<x + E l(J fiJia = A A/2. 
The c-orbital represents the semiconductor conduction band 
with dispersion e^, lying at an energy e Q above the hydrogenic 
1 s impurity bound state (the /-orbital), and V is the hybridiza- 
tion between them. Sj is the SU(A) spin operator of the fc- 
orbital. The hopping between the hydrogenic bound states [10] 
falls off exponentially with distance = jr, — rj\, = 
to exp(— Tij/a) for ry 3> a, the Bohr radius of the bound state. 
Consequently, the antiferromagnetic super-exchange coupling 

Jij = Jbexp(-2ry/a), (2) 

where Jo ~ (see footnote '). The projected Hilbert space of 
the /-orbital can be treated in the slave-boson formalism fj = 

bificr enslaved to a constraint on each site J2 a fi a fio + "i ^ = 
A/2. 

In this Letter, we focus on the insulating side of the un- 

1 /3 

compensated DS po < p c (p c a w 0.25 for Si:P) where 

— 1/3 

the average inter-site distance A = p ' a, which im- 
plies that tij — > 0. In this limit, the effective hybridization bV 
goes to zero as b — > 0, and the electrons become Mott local- 
ized on singly-occupied /-orbitals. This results in an effective 
Heisenberg Hamiltonian for the insulating uncompensated DS, 
^ = JW ~52i^j Jij^i ' Sj. The magnetic behavior of such a 
disordered Heisenberg system was largely explained by Bhatt 
and Lee via the SDRG method. Here we investigate the system 
within the large- A theory [8, 11], which leads to an effective 
mean-field Hamiltonian through the saddle-point approxima- 
tion, 

W = -^E^ (A£.A«+H.c.-|A«| 2 ), 

with the constraint J2 a floficr = N/2 (of self-conjugate spins) 
implemented through the local Lagrange multiplier Aj. Here, 

'it is well known (and expected) that corrections to exist due to 
anisotropy and other effects. For instance, in d = 3 there is an additional fac- 
tor of (rij/a) 5//2 multiplying Jy when ry S> a [10]. For our purposes, these 
corrections can be safely neglected in face of the highly disordered character of 
the dilute (insulating) regime. They only provide subleading (logarithmic) cor- 
rections as we confirmed numerically. Furthermore, this allows us to directly 
compare our results with those of the Bhatt-Lee theory [4], which also neglects 
them. 
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Fig. 1 : (Colour online) Normalized magnetic susceptibility of highly 
disordered 3D Heisenberg magnets evaluated with the Bhatt-Lee 
method [4], the large-iV self-consistent theory (for systems with No — 
512 spins (see footnote 2 ), and the geometric decimation procedure at 
concentrations n = ^-p a 3 = 0.004, 0.04, and 0.16. 

A-ij = 2J2a fiafj<?/N are valence bond (VB) operators and 

A jj = / Ay \ are variational parameters which minimize the 

free energy. They are solved self-consistently at A — > oo, for a 
given sample realization and temperature. The results are then 
averaged over 20 sample realizations. 

Numerical large-N solution. At any finite temperature, our 
large- A solution finds two types of spins: localized and VB 
spins. The localized spins are those isolated from all other ones, 
i.e., A = for all bonds connecting to them; their contribution 
to the magnetic susceptibility is simply a Curie term Xc(T) = 
fj/g/ksT. In contrast, each VB spin forms a singlet bond (A ^ 
0) with another spin; their contribution can be neglected at low 
temperatures. The low-T magnetic susceptibility is, therefore, 
well approximated by 

X (T) = p(T) Xc (T), (3) 

where p{T) is the density of localized/free spins at tempera- 
ture T. Figure 1 shows the normalized magnetic susceptibility 
x(T)/x(Jq) = Jop{T)/Tpo at concentrations n = ^-poa 3 = 
0.004, 0.04, and 0.16 (see footnote 2 ). The susceptibility di- 
verges at low temperatures, consistent with the SDRG results 
of BL [4]. This divergence is usually fitted by a power law in 
experiments, but we shall show later that it should be viewed as 
a logarithmic correction to the Curie law. The higher the doping 
concentration, the larger this correction since couplings among 
spins are stronger. At extremely low concentrations, all spins 
are essentially free and the magnetic susceptibility follows the 
Curie law. 

Geometric decimation procedure. The large- A ground state 
at zero temperature of such a highly disordered Heisenberg sys- 
tem is essentially a RS state, in which most spins form inert 
singlets (A = 1) with another spin and do not correlate with 
any other spin. To highlight this, we considered a simple four- 
spin cluster with antiferromagnetic couplings Jy > 0, and 

2 We have carefully verified that all our numerical results are robust with 
respect to finite-size effects. 
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J23 3> Jij for all 7^ (2,3). The large- N calculation 

shows that for T > J23, all bonds are zero and all four spins 
are free. As we lower the temperature to J23, spins S2 and S3 
start to form a VB singlet, A23 ^ 0, and no longer contribute 
to \ {T). Further reducing the temperature to J14, spins Si and 
S4 form another VB singlet. There is no resonance between 
the (2,3) and the (1,4) VB singlets. In contrast to the Bhatt-Lee 
SDRG method, in which there appears a renormalized cou- 
pling between SU(2) spins connected to a strong singlet pair, 
this effect can be shown to be of order 1 /N between SU(iV) 
spins [12], and thus drops out in the large-iV limit. While this 
simplification makes our large-TV model amenable to closed 
form solution, we shall demonstrate that it hardly affects the 
quantitative predictions of the model within the experimentally 
relevant temperature range (as shown in fig. 1). 

This also allows us to state a very simple geometric deci- 
mation procedure. We (i) search for the most strongly cou- 
pled spin pair, or equivalently, the shortest one (see eq. (2)), 
(ii) remove it from the system by coupling the spins in an in- 
ert singlet, and (iii) repeat steps (i) and (ii) until the desired 
energy (temperature) scale is reached. We should emphasize 
that no other renormalizations are involved during this deci- 
mation procedure. The density of free (undecimated) spins in 
eq. (3) is then given by p(T) = p a J T Q{J)dJ, where Q( J) is 
the distribution of the decimated couplings, shown in fig. 2(a) 
for ?i = 0.16. The distribution of nearest neighbor couplings, 
P{ J), is also plotted for comparison. Note the dramatic dif- 
ference between P( J) and Q(J) which stems from the fact 
that, during the decimation procedure, longer-distance nearest 
pairs are unavoidably generated. Therefore, Q( J) will always 
be singular yielding the divergence of x(^) at l° w tempera- 
tures. As depicted in fig. 1, this simple geometric decimation 
procedure captures the essential physics of the large-iV theory 
in describing the magnetic susceptibility of strongly disordered 
Heisenberg spin systems. 

Analytic solution. The geometric decimation procedure will 
give us a long-sought analytic description [13, 14] of the mag- 
netic properties of insulating DS if one can keep track of p 
as a function of the energy scale Q, = max{ Jy} (defined as 
the coupling to be decimated) or, equivalently, the length scale 
L = minjry •} (the distance between the spins in the pair to 
be decimated). Although the pair approximations [13] consid- 
erably simplify the calculations as compared to the SDRG and 
numerical cluster calculations, they fail to yield an analytic ex- 
pression for p. On the other hand, the analytic formula pro- 
posed by Ponomarev et al. involves a tunable parameter [14]. 
Here we present an accurate analytic solution without any ad- 
justable parameters for a general d-dimensional system. 

Since we remove hierarchically the closest spin pair, we can 
imagine each spin as a hard sphere of diameter L, which natu- 
rally incorporates the constraint that no spin pair is closer than 
L (see footnote 3 ). By removing the spheres that are touching 
each other, we continuously increase L until the next closest 
pair of spins touch each other. The rate equation governing the 



In the following, we neglect any other correlations beyond those imposed 
by the hard-sphere constraint. 




Fig. 2: (Colour online) (a) The distributions of the nearest neighbor 
couplings P(J) and of decimated couplings Q(J) at concentration 
n = 0.16. (b) Numerical results for the packing fraction 77 = pv as a 
function of decimation length L for d = 1,2, and 3; A = p 1 ^ d . (c) 
Comparison between the numerical and analytic (eq. (5)) results for 
the free spin density in the geometric decimation procedure for d — 3. 
Here, systems with No = 4 096 spins averaged over 3 000 samples 
were used (see footnote 2 ). 

density of free spins is given by 

dp = -2 d p 2 gdv, where g(p) = (1 - apv)(l - pv)~ d (4) 

is the radial distribution function [15] of a hard-sphere fluid: it 
gives the ratio of the density of particles at distance r by the 
mean density, given that there is a particle at the origin. Here, 
a is a constant which depends only on dimensionality (a = 0, 
0.436, and 0.5 for d = 1, 2 and 3, respectively) [15], and v 
is the excluded volume of each hard sphere. The negative sign 
comes from the fact that p decreases as L increases, and the 
decrease in p is proportional to the density of available spins p 
times the probability that two spins (hard spheres) touch each 
other, i.e., 2 d pgdv. The 2 d factor converts the radius of the hard 
sphere (raised to the power d) into its diameter. 

The solution of eq. (4) can be reduced to a quadrature, from 
which we can deduce that the packing fraction 7/ = pv in- 
creases monotonically with L, saturating at large length scales 
at r] c (~ 0.333, 0.182, and 0.0968, respectively, for d = 1, 
2 and 3). The results of a numerical solution of the decima- 
tion procedure are shown in fig. 2(b), from which we obtain 
77c ~ 0.2810(5), 0.156(1), and 0.082(2) for d = 1, 2 and 3 (see 
footnote 4 ). Since 77 <C 1 throughout the decimation procedure, 
our hard sphere liquid remains moderately correlated (away 
from the strong coupling regime in the vicinity to close pack- 
ing). This provides a dramatic simplification, since we are now 
well justified in using the virial expansion g _1 w 1 — (d — ajpv 
(this linearized expression is exact [15] in d = 1), and find a 



Small differences between the analytic and the numerical values of r/ c re- 
flect the higher order correlations we have neglected (see footnote 3 ). 



S. Zhou, J. A. Hoyos, V. Dobrosavljevic and E. Miranda 



closed form solution 

2 d 1P v = 1 - (p/po) 1 , with 7 = 1 + (d - a) /2 d , (5) 

which satisfies the initial condition p = po at v = 0. The 
magnetic susceptibility in eq. (3) is readily obtained by relating 
temperature and L via eq. (2), i.e., 2L = ahx(Jo/T). In the 
L, v — > oo (T — > 0) limit, the density decays asymptotically as 
p v^ 1 ~ Thus the magnetic susceptibility diverges at 
low temperatures according to 



X(T) 



Jo 



r[ln(J /T)] c 



(6) 



which can be viewed as a logarithmic correction to the Curie 
law instead of the power law divergence usually fitted to ex- 
periments. The free spin density p extracted from eq. (5) is 
plotted in fig. 2(c) as a function of L, in excellent agreement 
with the numerical result of the decimation procedure. There- 
fore, eq. (5) provides an accurate analytic solution, without any 
adjustable parameters, to the large-TV theory of the insulating 
DS. 

Comparison between SDRG and large-N. It is now natural to 
ask how reliable the large-TV theory is. To address this issue, we 
compare the well-known RS solution of the ID random Heisen- 
berg system obtained by the SDRG method [16, 17] with the 
analytic solution, eq. (5), of the large-iV theory. For randomly 
distributed spins, the length distribution of the nearest neighbor 
bonds is a Poissonian P(L) = po exp(-poL), which gives rise 
to a power-law initial coupling constant distribution 



Po(J) = 0(J)6(J - J) 



poa 
2Jo 



l-poa/2 



(7) 



In this case, the SDRG flow can be followed exactly through 
all energy scales, yielding [17] 



P = Po 



i + ^MJo/n) 



= Po (1 + PoL)~ 



(8) 



where the prime is added to distinguish this SDRG density from 
the large-TV result in eq. (5). In the asymptotic L — > oo limit, 
p' ~ L~ 2 , different from the L~ d behavior of the the large- N 
theory as shown in fig. 3(a). However, upon close inspection, 
the L dependences of p and p' (see fig. 3(b)) reveal that the 
breakdown occurs only above a length scale L* = 1/po = A, 
corresponding to a breakdown temperature 



T* = J cxp(-2A/a) 



(9) 



below which the renormalized couplings become important in 
the SDRG procedure. Above T*, however, the SDRG the- 
ory can be reduced to the simple geometric decimation proce- 
dure. The smaller the concentration p , the lower T* is. Since 
T* concerns only the energy scale at which the renormalized 
couplings become important, eq. (9) straightforwardly holds in 
higher dimensions. Interestingly, this result implies that a large 
class of highly disordered systems can be described by the ran- 
dom singlet picture above T* even though their ground states 
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Fig. 3: (Colour online) Comparisons between the densities of free (un- 
decimated) spins as functions of the length scale L obtained by the 
SDRG (p 1 ) and the large-A (p) methods, i.e., the geometric solution, 
in d = 1. The vertical dotted line highlights the breakdown tempera- 
ture T* (47 mK for ap c = a/A c = 0.25 (see footnote 5 )). 



are completely different (as they are in Refs. [18]). For in- 
stance, T* w 47 mK when ap 1 ^ = apl /3 = a/A c = 0.25, 
assuming Jo = 140 K from Ref. [19] 5 . Remarkably, the tem- 
perature window relevant for experiments is above the break- 
down temperature (left of the dotted line in fig. 3), which also 
explains the success of the BL theory. 

Finally, we would like to call attention to a caveat on eq. (6). 
As shown in fig. 3, the experiments take place in a tempera- 
ture range above T* in which both the SDRG and the geomet- 
ric decimation solutions coincide and before their asymptotic 
regimes have been reached. It is thus very clear that either the 
numerical solution or the analytic one in eq. (5) compare well 
with experiments. In d = 1 and above T*, eq. (5) can be well 
approximated by eq. (8). Again, it is clear that the apparent 
power-law divergence of the susceptibility seen in experiments 
should instead be interpreted as a logarithmic correction to the 
Curie law. 

Summary and outlook. We have shown how a variational 
large-A^ method provides a physically transparent and quanti- 
tatively accurate description of inter-site spin correlations on 
the insulating side of DS. In the presence of strong positional 
disorder, each localized spin forms a VB singlet with a rather 
uniquely defined partner, allowing for a closed-form solution 
of the problem in the large-A^ limit. 

Even more importantly, this approach opens a very attractive 
avenue to describe the behavior across the MIT. It is known that 
the large-A^ RVB approach correctly describes the high density 
Fermi Liquid state [8]. As we established that it also works in 
the opposite (insulating/Bhatt-Lee) limit, then it will also pro- 
vide a valid description of the transition by examining the two- 
orbital t-J model of eq. (1) with finite inter-site hopping ty. 
Each /-spin now has more than one choice: to still form a VB 



5 Corrections to Jij in eq. (2) may be important in order to compute the 
precise value of T*. If one naively inserts the factor of (rij /a) 5 / 2 , one gets 
that T* = J Q (A/a) 5 / 2 exp(-2A/a) instead of eq. (9), which increases T* 
by a factor of 32 at the critical concentration. However, note that A/a = 4 
is not much greater than 1. Thus, subleading corrections of order (r^j/a) 2 
[10] become important. To our knowledge, they are not known at the present 
moment. 
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singlet with another localized moment, or to undergo Kondo 
screening by conduction electrons. Similarly as in the large- N 
solution of the two-impurity Kondo problem [11] , we expect 
Kondo-screened sites to contribute to the formation of a coher- 
ent Fermi liquid, while VB singlet pairs to "drop out" from the 
conduction sea and remain Mott localized. Such gradual con- 
version of the correlated electron fluid into a localized VB solid 
may provide a microscopic underpinning for the phenomeno- 
logical "two-fluid" model [3] - possibly the key missing link 
for cracking the metal-insulator transition in doped semicon- 
ductors. 

We thank Ravin Bhatt for useful discussions, and Richard 
Stratt for very helpful comments about the hard sphere prob- 
lem. This work was supported by FAPESP grant 07/57630-5 
(EM), CNPq grant 305227/2007-6 (EM), and by NSF grants 
DMR-0506953 (JAH) and DMR-0542026 (VD). 

REFERENCES 

[1] Paalanen M. A. and Bhatt R. N., Physica B, 169 (1991) 
223. 

[2] SARACHIK M. P., in Metal-Insulator Transitions Revisited, 
edited by EDWARDS P. AND RAO C. N. R. (Taylor and Francis) 
1995. 

[3] Paalanen M. A. etal, Phys. Rev. Lett., 61 (1988) 597. 
[4] BHATT R. N. and LEE P. A., Phys. Rev. Lett., 48 (1982) 344. 
[5] MA S. K. et al, Phys. Rev.Lett., 43 (1979) 1434. 
[6] MOTT N. F., Metal-Insulator Transition (Taylor and Francis, 

London) 1990. 
[7] Anderson P. W., Science, 237 (1987) 1196. 
[8] AFFLECK I. and Marston J. B., Phys. Rev. B, 37 (1988) 3774; 

MARSTON J. B. and AFFLECK I., Phys. Rev. B, 39 (1989) 

11538. 

[9] For a recent review, see MIRANDA E. and DOBROSAVLJEVIC 

V., Rep. Prog. Phys., 68 (2005) 2337. 
[10] ANDRES K. et al.„ Phys. Rev. B, 24 (1981) 244; HERRING C. 

and FLICKER M., Phys. Rev., 134 (1964) A362. 
[1 1] JONES B. A. et al, Phys. Rev. B, 39 (1989) 3415. 
[12] HOYOS J. A. and MIRANDA E., Phys. Rev. B, 70 (2004) 

18040KR). 

[13] ROSSO M., Phys. Rev. Lett., 44 (1980) 1541; THOMAS P. A. 

and ROSSO M., Phys. Rev. B, 34 (1986) 7936. 
[14] PONOMAREV I. V. et al, Phys. Rev. B, 60 (1999) 5485. 
[15] SONG Y. et al, J. Phys. Chem., 93 (1989) 6916. 
[16] FISHER D. S., Phys. Rev. B, 50 (1994) 3799. 
[17] IGLOI F. etal, Phys. Rev. Lett., 86 (2001) 1343; IGLOI ¥.,Phys. 

Rev. B, 65 (2002) 064416; HOYOS J. A. et al, Phys. Rev. B, 76 

(2007) 174425. 

[18] LIN Y.-C. et al, Phys. Rev. B, 68 (2003) 024424; LAFLOREN- 

CIE N. etal, Phys. Rev. B, 73 (2006) 060403(R). 
[19] HlRSCH M. J. etal, Phys. Rev. Lett., 68 (1992) 1418. 



